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Abstract 

In this work, we consider complementary lattice arrays in order to enable a broader range of designs for coded aperture 
imaging systems. We provide a general framework and methods that generate richer and more flexible designs than existing ones. 
Besides this, we review and interpret the state-of-the-art uniformly redundant arrays (URA) designs, broaden the related concepts, 
and further propose some new design methods. 


Index Terms 

Coded aperture imaging, Golay complementary sequences, lattices, uniformly redundant arrays. 


I. Introduction 

I MAGING using high-energy radiation with spectrum ranging from V-ray to 7 -ray has found many applications including 
high energy astronomy HI, ||2l and medical imaging |l3l-||3. In these wavelengths, imaging using lenses is not possible 
since the rays cannot be refracted or reflected, and hence cannot be focused. An alternative technique to do imaging in this 
spectrum is to use pinhole cameras, in which the lenses are replaced with a tiny pinhole. The problem in these cameras is that 
the pinholes passes a low intensity light, while for imaging purposes, a much stronger light is needed. Increasing the size of 
the pinhole cannot solve this problem as it increases the intensity at the expense of decreased resolution of the image. Coded 
aperture imaging (CAI) is introduced to address this issue by increasing the number of the pinholes. A coded aperture is a 
grating or grid that casts a coded image on a plane of detectors by blocking and unblocking the light in a known pattern, and 
produces a higher signal to noise ratio (SNR) of the image while maintaining a high angular resolution ii, Q. The coded 
image is then correlated with a decoding array in order to reconstruct the original image. The deployment of pinholes and 
the decoding array are usually jointly designed such that perfect or near-perfect reconstruction is possible. Figure [T] gives the 
schematic diagram of a CAI system. 
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Fig. 1. An illustration of the CAI system 

Since the coded aperture is usually defined based on an integer lattice, it can be modeled as a two dimensional array. In 
general, we define the n-dimensional array C of complex-valued entries of size Li x • • • x 

L^[ci, * * * , Cn\ — 0, Vci, * * * , Cji G (ci, * * * , Cn) ^ [0, L\ 1] X • * • [0, Lyi 1]. 

For simplicity, C[ci, • • • , c„] is also denoted by C[a], where a = [ci, • • • , c„]^ G K". The decoding array D can be similarly 
defined. The set from which the elements of the aperture arrays take values from is referred to as an “alphabet”. In coded 
aperture imaging, a physically realizable coding aperture usually consists of 0/1 binary alphabet (representing closed/open 
pinholes). If multiple coded images are obtained with different aperture masks and the resulting digital projection images are 
suitably combined, a complex-valued array C becomes applicable 0, ||9l- For example, a {±1} coding aperture could be 
obtained computationally from two masks with openings at C’s 1,-1 locations. For an N-phase alphabet, it is calculated that 
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the number of {0,1} masks needed is {3N— l)/2 for odd N, and N for even N. The calculation follows from the observations 
that a root of unity in the form of a; + iy, xy 0 {i^ = —1) requires two masks, a pair in the form of (x + iy^x — iy), xy ^ 0 
require only three masks (a mask corresponding to x is shared), a pair in the form of {x + iy, —x + iy),xy 0 require only 
three masks (a mask corresponding to iy is shared), and l,—l,i,—i each requires one mask. Moreover, the development of 
hardware technology, e.g. spatial light modulators, may lead to realizable complex-valued physical masks. If both the coding 
and decoding systems use such masks, an analog reconstruction could be obtained. Due to the above reasons, we assume that 
the elements of an aperture could be unimodular complex numbers. 

For a planar object that is projected onto a gamma camera through the coding aperture, it can be shown that the object 
is perfectly decoded if C * D is a multiple m of the discrete delta function (5[r], where + denotes the convolution and (5[r] 
corresponds to an array with one centered at the origin and zero elsewhere 0, Q. The value of m is also called “the SNR 
gain”. Clearly, larger m is better, and m <uj, where uj = Li ■ ■ ■ Ln (the number of pinholes). Designing C and D is the key 
part of designing CAL 

It is useful to observe that the designs of apertures are intimately related to the concept of “autocorrelation”, and there 
are two typical types of “autocorrelation” for an array. Though the two types of autocorrelations are different, they can be 
both applied to the design of apertures through different approaches, as will be pointed out in this paper. One is “aperiodic 
autocorrelation”. The aperiodic autocorrelation function is given by 

{vi,- ■ ■ ,Vn) = ^ ,Cn]C'[ci,c„ + t;„], t;i, • • • ,u„ € Z, 

Cl,-" 

where c is the complex conjugate of c. The other is “periodic autocorrelation”, which is defined later in the paper. If we choose 
D = C~ , where C~[ci, ■ ■ ■ , c„] = C'[—Ci, • • • , —c„], then C * D gives the autocorrelation of C. 

Non-redundant arrays (NRA) have been introduced for arranging the pinholes in CAI, since they have the property that 
the aperiodic autocorrelations consist of a central spike with the side-lobes equal to one within certain lag (range of the 
argument vi, - ■ ■ ,Vn) and either zero or unity beyond the lag ifTol . Pseudo-noise arrays (PNA) ifTTll are another alternative, 
whose periodic autocorrelations consist of a central spike with —1 side-lobes, which lead to designs of a pair of arrays such 
that their convolution is a multiple of the discrete delta function mi. Twin primes, quadratic residues, and m-sequences are 
examples of PNA designs. NRA and PNA based designs are both referred to as uniformly redundant arrays (URA) 13, m, 
mi. However, the size of the URA structures is restricted and cannot be adapted to any particular detector El, m. Besides 
this, the SNR gain for URA is limited to a;/2 Q, ifTSl - lfTTl . Other designs that have also been used in CAI are geometric 
design ifTsl and pseudo-noise product design HU, but they are also available only for a limited number of sizes—the former 
design is for square arrays, and the latter one requires that pseudo-noise sequences exist for each dimension. 

Though it is generally hard to find a single pair of coding and decoding arrays, it may be easier to find several pairs which act 
perfectly while combined together. Based on this idea, we look for a broader range of designs of the coding arrays in this paper. 
We show that the aperture can then be customized to any shape on any lattice, meeting various demands in practical situations. 
Our work is inspired by Golay complementary arrays, which are defined as a pair of arrays whose aperiodic autocorrelations 
sum to zero in all out-of-phase positions. They have been used for pinhole arrangement in order to obtain the maximum 
achievable SNR gain, while eliminating the side lobes of the decoded image 0. We note that there is a natural mapping 
between a pair of Golay complementary arrays, say Ci and C2, and a CAI system consisting of two parallel coding/decoding 
apertures. It is illustrated in Figure |2] where Di = = C^. When an object goes through the system, the side lobes 


Cl 

C2 



Fig. 2. A CAI system with two parallel channels, with coding apertures Ci, C 2 , and decoding apertures C-^ , C 2 

are completely canceled out by the addition of the two decoded images. 

Though an aperture is usually defined based on an integer lattice, we consider the design problem in the context of a general 
lattice, naturally arising from practical implementations. For example, usually the distance between two pinholes should be no 
less than a given threshold due to physical constraints. It has been shown by L. F. Toth EOl that the lattice arrangement of 
circles with the highest density in the two-dimensional Euclidean space is the hexagonal packing arrangement, in which the 
centers of the circles are arranged in a hexagonal lattice. Thus, given the minimal distance allowed among pinholes, the most 
compact arrangement (thus with the largest possible SNR gain) is to arrange them on a hexagonal lattice. 

The outline of this paper is as follows. In Section |II] we briefly present related work on Golay complementary arrays (based 
on aperiodic autocorrelation), and then propose “complementary lattice arrays” and other related new concepts such as the 
“complementary array banks”, including Golay complementary array pairs as a special case. This general framework naturally 
leads to the new concept of “multi-channel CAI system” which extends the classical CAI system. We provide the concept. 






3 


theory, and the design framework. Due to the reasons mentioned before, our examples are based on two-dimensional hexagonal 
arrays and unimodular alphabets (which consist of unimodular complex numbers). Nevertheless, the methodology given in this 
work could be further generalized. In Section|III]we review the URA literature that is mostly based on periodic autocorrelations. 
We further generalize the related concepts in Section |IV] in a similar fashion. This leads to a new class of aperture designs, 
which have the desirable imaging characteristics of URAs, yet exist for sizes for which URAs do not exist. In Section |V] we 
provide computer simulations demonstrating the performance of our schemes. 


II. Concept, Theory and Design of Complementary Lattice Array 
A. Golay Complementary Arrays 

In this section we briefly introduce related works on Golay complementary array pairs. Golay lED first introduced Golay 
complementary sequence pairs in 1951 to address the optical problem of multislit spectrometry . Later, they were also 
used for many other applications, including horizontal modulation systems in communication ll 22 l . power control for multi¬ 
carrier wireless transmission ll2^ . Ising spin systems ll24l . channel-measurement ll25l . Il26l . and orthogonal frequency division 
multiplexing Ell. 

A concept that is related to complementary array pairs is the Barker array. It is a {±1} binary array C such that 

< 1, Vui,G Z,(z;i,^ (O,--- ,0). (1) 


Barker arrays are scarce. For n = 1, the known valid lengths are only 2, 3, 4, 5, 7, 11, and 13. For n = 2, it has been proved 

1 1 ’ 

1 -1 


that there is no Barker array for Li > 1 ,L 2 > 1 except 


. Another related concept is the NRA, which also satisfies 


the condition ([T]). Its only difference compared with Barker array is that it is {0, l}-binary. 

Golay complementary array pairs address the scarcity of Barker arrays and NRAs. The basic idea of Golay complementary 
array pairs is to use the nonzero part of one autocorrelation to “compensate” the nonzero counterpart of the other ED. 
Specifically, a pair of arrays Ci and C 2 of size Li x • • • x is a Golay complementary array pair, if the sum of their 
aperiodic autocorrelations is a multiple of the discrete delta function, i.e. 


,Vn) + - ■ ■ ,Vn) = 0, 


Vui, • • • , Un G Z, (Ui, • • • , Un) ^ (0, • • • , 0). 


The initial study of Golay complementary sequence pairs (n = 1) was for the binary case. Binary Golay complementary 
sequence pairs are known for lengths 2 , 10 ESj and 26 ESl . It has been shown that infinitely many lengths could be synthesized 
using the existing solutions ll2^ . Specifically, binary Golay complementary sequence pairs with length 2^G0^^26^^ exist, where 
ki,k 2 , fcs are any nonnegative integers. Besides, no sequences of other lengths have been found. Later on, larger alphabets were 
considered, including 2”-phase EOl. A-phase for even N |[3T1 , the ternary case 2 l = {— 1 , 0 , 1 } |[32l-|[34l, and the unimodular 
case ES. Here, an alphabet 21 is called A-phase, if it consists of Ath roots of unity, i.e. 21 = {^ : = 1}; it is unimodular 

if 21 = {C : Id = !}■ 

In 1978, Ohyama et al. El constructed binary Golay complementary array pairs (n = 2) of size 2^^ x 2^^. The size is then 
generalized to 2''G0^^26^^ x where kj,j = 1 , ■ ■ ■ ,6 are any nonnegative integers E3, ED. 

We look for broader concepts and designs than complementary array pairs. The examples provided in this paper are for the 
two-dimensional case, but they can be easily generalized to higher dimensions. We start with the definitions in the following 
section. 


B. Definitions and Notations 

Definition 1. A lattice in K” is a subgroup o/R” which is generated from a basis by forming all linear combinations with 
integer coefficients. In other words, a lattice £ in R" has the form 

£= gZ 

U=i 

where {edidi forms a basis o/R". 

For example, the integer lattice 1? is generated from the basis ei = (1, 0), 62 = (0,1). The hexagonal (honeycomb) lattice 
A 2 is generated from the basis ei = ( 1 , 0 ), 62 = (— 5 , ^)- 

A classical array is based on an integer lattice. We now give the definition of an array that is based on a general lattice. 

Definition 2. Let 2 be a lattice. A lattice array A defined over alphabet 21 and with support LI is a mapping C[-\ : £ —?► 21, 
such that C[a] = 0 for all a LI and C[a\ G 'Qi for all a £ LI. The number of the elements of Ll (array size) is denoted by 
|0|. Bb denote hy (J when there is no ambiguity. In other words, 

C[a]=0, yaG2\Ll, 
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where C[a\ is the entry at location a. 

The following terms are made to simplify the notations. 

• Define as the shifted copy of by t (for t G £), if 

= C'^’^’®[a-t],Va G £. 

For brevity, {f\ is simplified as C{f\. 

• Assume that two arrays based on the same lattice £, but not necessarily on the same area. 

The addition of Ci and Ci, C = C\ is an array whose entries are the addition of corresponding entries in C\ and 

Cl, i.e. 

f2 = f7iUf22, C[a] = Ci[a] +C 2 [a\, \/a G fl. 

• A set of arrays non-overlapping if 

n Vtm 2 =0, Vmi, 7712 G 1, 2, • • • , M, mi ^ m 2 . 

Definition 3. Assume that the lattice £ is generated from {ei}2-i. The aperiodic autocorrelation function is 

,i;„) = ^ C'[a]C[a + i;iei H-hWne^], vi,--- G Z. 

The aperiodic crosscorrelation function of two arrays Ci and C 2 is 

■ ■ ,Vn) = ^ Ci[a]C2[a + i’iei H-h z;„e„], vi,--- ,ii„ G Z. 

oGn 

Sometimes, and are respectively denoted by C * C~ and Ci * Cf. 

Definition 4. A complementary array bank consists of pairs of the crosscor¬ 

relations is a multiple of the discrete delta function: 

M M 

Y^Cm*Df,= Y, = u:S[r], 

m—1 m—1 

where ui is a constant. M is called the order, or number of channels. 

Remark 1. There is a natural mapping between a complementary array bank, say {{Cm, Dm)}m=i’ ^ CAI system 
consisting of M parallel channels, each of which consists of a pair of coding and decoding apertures (Fig.\^. When a source 
image comes, it is coded and decoded through M channels simultaneously, and is then retrieved by simply adding the decoded 
images from all the channels. The multi-channel CAI system proposed here provides a generalized solution to CAI design, by 
including a classical CAI system as a special case. In the remaining part of Section we mainly study the complementary 
array sets, which may provide insight into the theory and design of complementary array banks in general. 


Cl 


or 

C2 


( 

CM 
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Fig. 3. A CAI system with M paiallel channels, with coding apertures Ci, ■ ■ ■ , Cm, and decoding apertures 

Definition 5. A set of arrays is a complementary array set, if the sum of their aperiodic autocorrelations is 

a multiple of the discrete delta function, i.e. 

M 

A'^^{vi, ■■■ ,V„) = 0, yVi,V2 G Z, (ill, • • • ^ (0, • • • ,0). 

m—1 

A Golay complementary array pair is the special case when M = 2. 


( 2 ) 
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Remark 2. In practice, the pinholes on an aperture may only change the phase of a source point. Therefore, we assume a 
unimodular alphabet by default. It is clear that if a set of non-overlapping arrays are based on unimodular/N-phase alphabets, 
the addition of them is also based on an unimodular/N -phase alphabet. 

The autocorrelation of any array is the same as that of its shifted copy. This is because for any Vi, - ■ ■ , £ Z, t £ £, we 

have 

■ ■ ,Vn) = ^ C{t}[a\C{t}[a-\-viei-\ - \-Vne„] = ^ C[a - t]C[a - t-\-viei-\ -h 

= ^ C[a\C[a + viBi H-h , Vn). 

Furthermore, if {Cm}m=i complementary array set, then {Cm{tm}}m=i j £ ■C also forms a complementary array 
set. In other words, a complementary array set is “invariant” under shift operation. 

Based on a unimodular alphabet, a complementary array set {Cm}m=i satisfies 

M 

A^’^io,--- , 0 ) =M\n\. 

m—1 

Thus, the sum of the autocorrelations may be written as a multiple of the discrete delta function: 

M M 

YCm*C- = Y = AI\mr]- 

m—1 m—1 


C. Inspirations from Ohyama et al.’s Design 

In Ohyama et al.’s design, £ is an integer lattice, and the number of complementary arrays is M = 2. The design consists 
of two steps; 

• First, choose the following complementary sequence pair: 

Ci = [l 1], C'2 = [l -1]. (3) 

• Second, design complementary array pairs of larger sizes in an inductive manner. Assume that we already have a 
complementary pair Ci, C 2 , with Ci * Cf -\- C 2 * Cf = 2u!S[r], where w is constant. Let 

Cl = Ci{fi} + C2{t2}, C 2 = Ci{fi} - C2{f2}, (4) 


where the shifts ti and t 2 are arbitrarily chosen. 

The validity of the construction (|4]i is clear from the fact that 

Cl * Cf =Ci{fi} * Ci{fi}- + C2{f2} * C2{t2}- + Cilfi} * C2{f2}- + C2{t2} * Ci{fi}-, 
C 2 * df =Ci{fi} * Ci{fi}- + C2{f2} * C2{t2}- - Ciih] * C2{t2}- - C2{t2} * Ci{fi}-, 

and thus 


Cl + Cl + C 2 * C 2 


2 (Ci{fi} * Ci{fi}- + C2{f2} * C2{t2}~) = 2 (Ci * Cf + C2 * Cf) = 4 ujS[r]. 


(5) 


In practice, ti and <2 are chosen properly so that Ci{fi} and C 2 {f 2 } do not overlap, which guarantees that Ci and C 2 are still 
based on unimodular alphabets. For example, after applying Equation (|4]i to Q once, we have two possible complementary 
array pairs: 


-1 1 ]; ( 6 ) 


(7) 

The process of design is also illustrated in Fig. |4] By continuous application of the above design process, complementary pairs 
of size 2^1 X 2^'^ (for any nonnegative integers ki, ^ 2 ) can be designed. 

Inspired by the above design for complementary array pairs on integer lattices, we look for a “seed” (similar to (O) and 
a related scheme to “grow” the seed (similar to (IHi) for the design of complementary hexagonal arrays. Admittedly, we may 
build a simple mapping between two-dimensional arrays on a square lattice and a hexagonal lattice (or other lattices) below: 


Ci = [l 1 1 -1], C2 = [l 1 


or 


Cl = 


■ 1 1 ■ 


1 -1 

_1 

II 

(N 

^3 


1 ] 

-1 ] 


fjS. ,n .21 + C2 * 63] = C^ [ci * ej* + C2 * 62], Vci, C2 £ Z, 


( 8 ) 
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The first way: 



+ 

— 



- 


The second way: 



Fig. 4. Illustration of the design of complementary array pairs on integer lattices 


where the superscripts s and h respectively denote square and hexagonal lattices. Under the above mapping, a set of comple¬ 
mentary arrays on a square lattice are still complementary on a hexagonal lattice. This is due to the fact that the autocorrelation 
of an array is only with respect to the coefficients ci,C 2 . Nevertheless, the lattice array naturally arises from practical designs. 
Consider the scenario where a two (or three)-dimensional coded aperture is to be built that has pinholes arranged on a certain 
(suitably chosen) type of lattices which adapts to a particular physical aperture mask. The designs may preferably be based 
directly on that lattice instead of mapping to a square lattice (with zero elements padded in various regions) first and then 
mapping back to the original one. 

D. Design for The Basic Hexagonal Array of 7 Points 

We first study a very simple hexagonal pattern that may act as a “seed”. It is an hexagonal array of 7 points, which is shown 
in Fig. 12 After that, we consider possible ways to “grow” the seed. 



Fig. 5. Basic hexagonal array with 7 points 

We start from considering the existence of hexagonal complementary array pairs, i.e. the order M is 2. 

Theorem 1. For the basic 7-points hexagonal array, there exists no complementary array pair with unimodular alphabet (Fig. 

121 ) 



Fig. 6. There exists no such complementary pair with unimodular alphabet. 


The proof is given in Appendix One may be further interested in the existence of a hexagonal complementary array pair 
if the basic array does not have the origin 0 (Fig. |7]i. In fact, it does not exist, either. 

Theorem 2. For the array in Fig. [2 there exists no hexagonal complementary array pair with unimodular alphabet. 

The proof is given in Appendix |B] The non-existence of complementary array pairs for the array in Fig. 0 motivates us 
to further consider higher order M. We use the notation of “design parameter” for brevity. For a particular array pattern, if 
there is a complementary array set with M arrays and an A^-phase alphabet, the pair (M, N) is called its design parameters. 


































Fig. 7. Basic hexagonal array with 6 points 



Furthermore, if the array sizes are equal to L, we may refer to the triplet (M, iV, L) as its design parameters whenever there 
is no ambiguity. 

Fortunately, complementary array triplets with unimodular alphabet exist. In fact, we have found more than one designs 
with (M, N, L) = (3, 3, 7). The following is an example. 

Design 1. Let Q = exp(i 27 r/ 3 ). Let Ci = {xk}l^Q,C 2 = {yk}k=o^^3 = {^k}k=o denote the entries of three hexagonal 
arrays shown in Fig. |5] Then 


{^4Lo = {C^ C°, C^ C^ C°, C^ C°}, telLo = {c\ C°, C^ C\ C°, c'}, Uk}l=o = {c\ C°, C\ C\ C^ C°, c'} 


form a complementary array set (Fig. 0. 



Fig. 8. A complementai'y triplet with 3-phase alphabet, i.e. (M, N, L) = (3, 3, 7), where is represented by k , for fc = 0,1, 2 

We have also found more than one designs with (M, L) = (4, 2, 7). The following is an example. 

Design 2. Let Ci = {iCfc}®^Q,C '2 = {yk}k=o^^3 — {^k}k=o^^i = denote the entries of four hexagonal arrays 

shown in Fig. |5] Then 

{xk}Lo = { 1 , 1 , - 1 , 1 , 1 , - 1 , 1 }, {yk}l=o = {- 1 , 1 , 1 , - 1 , - 1 , 1 , 1 }, 

{^4Lo = {1>1,-1,1,-1,1,1}, {t«fc}Lo = {1,1,1,-1,1,-1,!} 

fonn a complementary array set (Fig. 0- 



Fig. 9. A complementary quadruplet with 2-phase (binary) alphabet, i.e. {M,N,L) = (4,2,7), where ±1 is represented by di 


E. Methodology for Designing Larger Arrays 

We now consider how to “grow” the seed that we have found in order to design more and larger arrays. If + and * are 
respectively considered as addition and multiplication operations, then (|4]i and (|5]l can be written in symbolic expression: 


where 


C = HC{t}, 

= C{f}'^H'^HC{f}- = 2C{f}^C{f}-, 


1-1 

, c- = 

' Cf ■ 
. ^2” . 

, H = 

■ 1 1 ■ 
1 -1 

, C{t} = 

■ Ci{F} ' 

C2{t2} 


(9) 

( 10 ) 


and is the conjugate transpose of H. The key that C remains to be a complementary pair is that H satisfies H^H = 21, 
i.e. iT is a Hadamard matrix. This observation could be generalized to the following result. 


Theorem 3. Let U = \ujnk\MxM be a unitary matrix up to a constant, i.e. U'^U = cl, where c > 0. Assume 
is a complementary array set. Then, {Cm}m=i is also a complementary array set, where 


M 

m = 1,2, • • • , M, 

fe=i 

■ ■ ■ tM are arbitrarily chosen, and u ■ C is an array that multiplies each entry of C by the scalar u. 
Proof. Define a vector space on sQ with the addition operation + defined in Definition |2] and the variables 



■ Cl ■ 


Ci{ti\ 

c = 


,C{t} = 



. Cm . 


Cui^M} 


Define a quadratic form with the multiplication operation * defined as the convolution. The sum of aperiodic autocorrelations 

of is 


M 

'^Cm*C-= C^C- = C{f}'^U'^UC{t}- = cC{f\'^C{t}- 

m—1 


M M 

— C ^ ^ — C ^ ^ Cm * — 

m—1 m—1 


cMuj5[r\. 


( 11 ) 

□ 


Remark 3. We may be more interested in the case where 

1 ) 21 is an N-phase alphabet; 

2) U is the Lourier matrix: Fm = [fmk]MxM, fmk = exp{i27r(m - l)(fc - 1)/M}; 

3) The shifted arrays C^’ fc = 1, • • • , M do not overlap. The complementary array set with design pa¬ 
rameters (M, |n|) fi/irifcl = |n|,A: = l,-- - ,M) becomes {Cm}m=i with design parameters {M, lcm(A^, M), M|n|) 

according to Theorem 3, where 1cm stands for least common multiple. Besides this, we have c = M, uj = |f2| in ( 1771) . 

Theorem [3 provides a powerful tool to design more complex complementary arrays, which will be illustrated in Subsection 
III-FI For future reference, we also include the following fact. 

Remark 4. Assume that • • • , are K complementary array sets on the same lattice, and the set has design 
parameters {Mk, Nk), V fc = 1, • • • ,K. Then, it is clear that C = U • • • U can be thought as a new complementary 
set with design parameters ^k, lcm(A'^i, • • • , Nk)). 

F. An Example—Design for the hexagonal array of 18 points 

The basic hexagonal array of 7 points studied in Subsection III-DI contains two layers, with 1 and 6 points, respectively. We 
now study the hexagonal array that contains one more layer (shown in Fig. fTOl) . The design for this hexagonal array is not 
very obvious, so we delete the center element §, i.e. layer 1. The remaining 18 points are grouped into six basic triangular 
arrays: Ci,C 2 , C 3 , Ci,C 2 , C 3 , which is shown in Fig. [TT] 

* It is not necessarily a field. 

“In fact, designs always exist for an arbitrarily shaped array, as will be discussed later. We delete the center point primarily for a cute solution. 
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Fig. 10. Hexagonal an'ay with three layers (19 points) 



This motivates the design of complementary array triplets for the basic triangular array shown in Fig. [12] where is 

represented by k , for k = 0,1,2, = exp(i27r/3). In Fig. [T2| {Ci, (72, Ca} form a complementary array set. Due to 
symmetry, its rotated copy, ( 7 a,( 72 ,( 73 , also form a complementary array set. Applying Theorem |3| we obtain the following 
design. 



Fig. 12. Basic triangular complementary array triplets with 3-phase alphabet, i.e. (M,N,L) = (3,3,3) 



Design 3. Let 


Cl 

Cj 



4 
. ^6 


= -F 6 


■ c[{ti} - 

Cl{t2} 

C2{t3} 

(72 {* 4 } ’ 

C'siU} 

. ( 73 {* 4 } . 
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Fig. 13. Complementary array set with parameter {M,N,L) = (6,6,18) 


where ti, • • ■ , such that Ci, • • • ,Cq are arranged to form a hexagonal array of 18 points, as is shown in Fig. [77] 

Then {Cm}m=i forms a complementary array set. The design is also shown in Fig. [75] where Cf is represented by k, for 
k = 0, - ■ ■ , 5, C = exp(z27r/6). 

G. Complementary Array Bank 

Up to this point, we have assumed that the coding array C and decoding array D are related via D\r\ = C[—r\. The design 
of CAI thus reduces to the design of complementary array sets. Then we extend the autocorrelation to crosscorrelation, and 
the design of complementary array sets is accordingly extended to that of complementary array banks. The following result is 
a generalization of Theorem |2 and its proof is similar to that of Theorem |3| 

Theorem 4. Let © = = {(j)rnk\ € be two matrices satisfying ©^$ = cl for some positive constant a For 

a given lattice £, assume {{Ck, is a complementary array bank. Then, {(Cm/Dm)}m=i olso a complementary 

array bank, where 

M M 

Cm — ^ ^ ^mk * Ckflk ^, Dm — ^ ^ fmk * Dk\tk\ , 771 = 1, 2, • • * , Af, 

k=l k=l 

- ■ ■ tM are arbitrarily chosen, and u ■ C is an array that multiplies each entry of C by the scalar u. 

Remark 5. We may be more interested in the case where 

1 ) have N-phase alphabets; 

2) © is equal to $ and it contains M orthogonal columns of the complex Fourier matrix Fj^ (thus M < M); 

3) The shifted arrays Ck{tk} do not overlap, neither do fc = 1, • • • , M. 

Furthermore, if we assume that Ck = Dj^, k = 1, • • • , M and the array sizes are equal to L, then the complementary array 
set has design parameters {M, N, L), and {Cm}m=i design parameters {M,lcm{N,M),ML). 

Lemma 1. Any single point, as the simplest array on any lattice, forms a complementary set. 

Remark 6. Lemma\I\is a trivial but useful result. It follows from Definition\^and the fact that the aperiodic autocorrelation 
of a single point is always the discrete delta function. By employing Lemma [7] Theorem [i] and Theorem [7] it is possible to 
design complementary arrays of various support LI. 
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Corollary 1. For an arbitrary set fl on a lattice, there exists at least one complementary array set with support for any 
order M such that M > |0|. 

Remark 7. [Augmentation using Theorem ^ Due to Theorem 0 we may let M > M for practical purposes. For example, 
we may choose M = 2"* ( for a positive integer m) such that that U is a ±1 Hadamard matrix and the growth of alphabet 
(N) could be well controlled. We call this “augmentation” procedure. Augmentation is important, because it is often desirable 
to reduce the size of the alphabet, and thus the cost of practical implementations. The following design is an example of 
augmentation. 

Design 4. We use several basic arrays to compose a smile face shown in Fig. [T^ The colors indicate different basic 
complementary array sets: The two green arrays (at the lower and upper boundaries) fonn a basic array set with parameters 
(M,N,L) = (2,2,8). So do the brown ones (at the lower left and upper right boundaries) and cyan ones (at the lower 
right and upper left boundaries). The two blue arrays (at the lower and upper boundaries) fonn a basic array set with 
(M,N,L) = (2,2,4). The two yellow arrays (single points at the lower and upper boundaries) fonn a basic array set with 
(M, N, L) = (2,1,1). The three black arrays (the eyes and nose) form a basic array set with (M, N, L) = (3, 3, 3) (Fig. I72I) .' 
so do the three red ones (part of the mouth). The four purple arrays (the rest part of the mouth) form a basic array set with 
parameters [M, N, L) = (4, 2, 3), which could be obtained by applying Theorem^with Remark\^and M = 4 to a set of three 
single-point arrays. By applying Remark^and Theorem\^to the 20 arrays, a smile design with {M,N,L) = (20,60,88) 
could be obtained. Due to Remark^ another smile design with {M,N,L) = (32,6,88) could be obtained. 



Fig. 14. Smile design, with (M, N, L) = (20, 60, 88) (without augmentation) or (32, 6, 88) (with augmentation) 


//. Design for Infinitely Large Hexagonal Arrays 

By choosing a proper seed and growth scheme, we may be able to design infinitely large hexagonal arrays. The following 
is an example. We first design a complementary array set with M = 7, as a seed. 

Design 5. The union of Design\I\with design parameters {M,N,L) = (4,2,7) and Design^with (M,N.,L) = (3,3,7) is a 
design with (M, N, T) = (4 + 3, lcm(2, 3), 7) = (7, 6, 7) (Fig. 13), based on Remark^ 

Design 6 . By repeated application of Theorem^with U = Fz to Design^ we obtain a design with parameters {M,N,L) = 
(7,42, 7^) for any positive integer L The design is illustrated in Fig. [75] where the colors indicate the process of “growth”. 
In fact, applying Theorem\^to Design^ once (with Remark\^ conditions), we obtain a larger array set with (M,N,L) = 
(7, lcm(6, 7), 7x7) = (7, 42, 7^). Fig. \15(a)\ illustrates how the 7 arrays of size 7 (indicated by different colors) are combined 
to form larger arrays. Similarly, applying Theorem\^to the {M,N,L) = (7,42,7^) design once, we obtain a larger array set 
with (M, TV, L) = (7, lcm(42, 7), 7 x 7^) = (7,42, 7^). Fig. \15(b)\ illustrates how the 7 arrays of size 7^ (indicated by different 
colors) are combined to form larger arrays. Further applications of Theorem\^will not increase M,N, but will increase L. 

As an alternative, the following design is also for 7^-point hexagonal array, but with different elements. 

Design 7. Bfe keep applying Theorem |i| with U = Fz to a single-point array, e.g. with entry 1, we obtain a design with 
parameters (M, N, L) = (7, 7, 7^) for any positive integer L To see how it works, first consider a complementary array set 
with (M, N, L) = (1, 1, 1) (a single-point array). Taking the union of 7 such array sets as in Remark^leads to an array set with 
(M, N, L) = (7,1,1). Then, applying Theorem\^once leads to an array set with (M, TV, L) = (7, lcm(7,1), 7x1) = (7, 7, 7). 
Further applications of Theorem |3 will increase L, but not TVf, TV. The design could also be illustrated by Fig. [73 
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(b) Illustration of one of the 7 arrays of size 7^ obtained as a result of applying 
Theorem O to a complementary array set of 7 arrays of size 7^ (indicated by 
different colors) once 


(a) Illustration of one of the 7 
arrays of size 7^ obtained as a 
result of applying Theorem [3] to 
a complementary array set of 7 
arrays of size 7 (indicated by dif¬ 
ferent colors) once 


Fig. 15. Complementary array set with parameter (M, N) = (7, 42, 7^) (Design[6j, or (M, N) = (7, 7, 7^) (Design!?) for any positive integer -f 


III. URA, HURA AND MURA Based on Periodic Autocorrelation 
I n this section, we review related works on URA, including HURA and MURA. 

A. URA 

We first introduce the concept of periodic autocorrelation and “pseudo-noise” that are important to the design of URA. 
Definition 6. Let C = {C[ii, • • • ,fn]} be an infinite array on an integer lattice, which satisfies 

,in] = Vii,--- e = Ji (mod Li),--- = /„ (mod L„), 

Ll, - ■ ■ , Ln, /l, • • • , G N, (Ji, • • • , In) G [0, Li — 1] X • • • X [0, Ln — 1]. 

The finite array within [0, Li — 1] x • • • x [0, — 1], denoted by c, is called the basic array. C is called the periodic extension 

of s. The periodic autocorrelation function of C (or c) is 

,Vn) = ^ ,in]C[il+Vi,--- ,in + Vn], Vi,---,VnGZ. 

iltz[0,Li — 1 ],‘" — 1 ] 

The periodic crosscorrelation between two arrays are similarly defined. 

A section of an infinite array C would be a valid URA aperture, if there exists a finite array D such that C * D~ is a 
periodic extension of the discrete delta function. For a detailed discussion about the benefits and implementations of periodic 
extension, please refer to 121 . 

Definition 7. An array of size Li x ■ ■ ■ x Ln is a pseudo-noise (PN) array if 

(1) it is {zLl}-binary; 

(2) all out-of-phase correlations are —1 MIV . i.e. 

A^{vi,- ■ ■ ,Vn) = -1, Vi,-- ■ e Z, Ui ^ 0 (mod Li), • • • , or Vn ^ 0 (mod Ln). 
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In 1967, Calabro and Wolf IfTTI showed that a class of two-dimensional PN arrays could be synthesized from quadratic 
residues. The arrays are of size pi x p 2 , where pi,p 2 are any prime numbers satisfying p 2 — pi = 2 


D[ii,i2] = 

where {i/p),i G Z is Legendre operator: 

{i/p) = 


— 1 *2 = 0 mod p 2 

1 *1=0 mod pi, *2 ^ 0 mod p2 

(*i/pi)(* 2 /p 2 ) otherwise 


0 

1 

-1 


* = 0 mod p 

3x^0 mod p, s.t. i = mod p 
otherwise. 


In 1978, following from the above result. Cannon and Fenimore Q designed C and D such that C * D is a periodic 
extension of the discrete delta function. The design is given below, where pi ,p 2 is a twin prime pair. The coding array is C: 

1 (*iM)(*2/P2) = 1 

0 *2 = 0 mod p 2 

1 *1 = 0 mod pi, *2 ^ 0 mod p 2 

0 otherwise. 


C'[*l,*2] = < 


( 12 ) 


The decoding array of C is 79 , where D is: 

79[*i,*2] = 


1 

-1 


if C'[*i,* 2 ] = 1 
if C'[*i,* 2 ] = 0 . 


It is shown that 


C,D- = ™^ 


(5[r] (within one period). 


(13) 


URA may also be designed from Maximal-length Shift-register Sequences or m-sequences 1381 . m-sequence is another class 
of PN sequences, which have lengths n = 2^ — 1 with k being any positive integer. They are sometimes referred to as “PN 
sequences” or “m-sequences” 13^ . In 1976, MacWilliams and Sloane showed how to obtain PN arrays from m-sequences 
l39l . Let S be an m-sequence of length n = 2^ - 1. If n = 7*17*2 such that tii and 7*2 are relatively prime, a PN array H is 
designed below: 


^^[*11*2] = S[i], where * = *1 mod 7*1, 0 < *1 < **1, and * = *2 mod 7*2, 0 < *2 < 7*2. ( 14 ) 

We note that when sum( 7 T) = — 1 ( l39l . Property IP-/L*), we can design a URA with coding array C = (—77 + J)/2 and 
decoding array 79 = —H~ on integer lattices, where J is a unit array, i.e. with all elements equal to one. 


B. HURA 

In 1985, Finger and Prince na designed a class of linear URA, i.e. sequences C and 79 such that C * 79 is a periodic 
extension of the discrete delta function. The design is based on PN sequences which in turn come from quadratic residues. 
Then, by mapping linear sequences onto hexagonal lattice, they proposed the hexagonal uniformly redundant arrays (HURA). 
In the first step, they constructed the following sequence of length p, where p = 3 mod 4 is a prime: 

1 if * = 0 

— {i/p) otherwise 


79[*] = 


Let C = {D + J)/2. Then the following identity holds: 

C*D- = 


(5[r] (within one period). 


In the second step, they map the sequence 79 onto a hexagonal lattice: 

77'[*iei -f * 262 ] = 79[*i -f t* 2 ], 


(15) 


(16) 


where r is an integer to be chosen. H is called the Skew-Hadamard URA. It is easy to see that the correlation between H 
and (77 + J)/2 is a multiple of the discrete delta function, just like the one-dimensional case. 

As to the choice of lattice and r, it is well stated in CSl that “The freedom available in this procedure rests in the choice of 
the lattice, the choice of the order p, and the choice of the multiplier r. The lattice type will determine what symmetries can 
occur • • • The multiplier r determines the periods of the URA and hence the shape of the basic pattern.” Furthermore, HURA 
are those with hexagonal basic patterns, when the lattice is chosen to be hexagonal. The qualified p is either 3 or primes of 
the form 12 fc + 1 ifTSl . 

Besides the fact that HURA are based on hexagonal lattices, they are antisymmetric upon 60 degree rotation. This property 
provides for effective reduction of background noise El, G). Due to similar reasoning, the designs proposed in Section HI] also 
obtain robustness against background noise. 
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C. MURA 


It has been shown that PN sequences, together with the URA and HURA that are based on them, could be made with prime 
lengths of the form 4fc + 3. Gottesman and Fenimore M proposed the modified uniformly redundant arrays (MURA), which 
further increased the available patterns for CAI. MURA exist in lengths p = 4fc + 1 where p is a prime. 

The design of MURA also starts with a sequence D which is then mapped onto a hexagonal lattice, following the same 
procedure as HURA. Recalling URA and HURA designs from Subsections IIII- Al and IIII-Bl we may design using the procedure 
below: 

Step 1. Let D he a PN sequence (array); 

Step 2. Let the coding array C he {D + J)/2, and the decoding array be D~; 

Step 3. (optional) We map sequences onto a two-dimensional lattice (see Equations (fT4l) and (fThl) ). 

However, the design of MURAs is less straightforward, because D is not a PN sequence and C ^ {D + T)/2. One way to 
design MURA sequences is: 


CV\ 


m 


0 j = 0 mod p 

1 3x^0 mod p, s.t. i = x^ mod p 

0 otherwise 

1 z = 0 mod p 

1 C[i] = 1, z ^ 0 mod p 
— 1 otherwise. 


It is easy to verify that for any z; ^ 0 mod p, we have J2i=o C![i]D[i + u] = 0. 

Gottesman and Fenimore also gave a class of MURA for integer lattices. The coding array is the same as (HU, except for 
a change of the size: pi = P 2 = p. The decoding array is D~, where 


D[ii,i2] 


1 zi + Z 2 = 0 mod p 

1 C'[zi,i 2 ] = 1 ,U + *2 = 0 mod p 

— 1 otherwise. 


IV. New URA Constructions 
A. URA from Periodic Complementary Sequence Set 

In this section, we first briefly summarize some similarities and differences between the aperiodic-based and periodic-based 
designs of CAI, and then propose a new design framework that is based on periodic autocorrelation. 

In the aperiodic case, the elements of arrays are assumed to extend only over some finite area and be zero outside that area. 
This fact provides great convenience for the design of complementary array sets/banks, since several arrays could be easily 
concatenated while maintaining the unimodular alphabet during the “growth” process. In addition, the concept of “bank” and a 
growth scheme make the aperiodic-based designs more flexible. For example, we have shown how to make CAI aperture with 
arbitrary patterns. In the periodic case, the arrays were assumed to be periodic and infinite in extent. The resulting correlations 
are calculated over a full period. The usual way to design is to first design sequences with good autocorrelation property, 
e.g. pseudo-noise, and then map them onto arrays. As to practical implementations, periodic-based designs often require the 
physical coding aperture to be periodic extensions of some basic patterns to mimic the periodicity, while aperiodic-based ones 
do not. 

Despite their differences in principles and implementations, the idea of “complementary” can also be associated with periodic 
correlations, leading to the following concept that is similar to complementary array sets in Section HI] 

Definition 8. A set of arrays with the same basic pattern is a periodic complementary array set (PCAS), if the sum of their 
periodic autocorrelations is a periodic extension of the discrete delta function. A one-dimensional PCAS is also referred to as 
a periodic complementary sequence set (PCSS). The notation “design parameters” {M,N) or {M,N,L) is similarly defined 
as in Subsection 1//-D1 

As discussed before, URA (including HURA) require the lengths of sequences to be prime numbers or 2^ — 1, so the possible 
sizes of URA arrays are quite limited. However, the above concept produces more admissible lengths, offering more choices 
in selecting an aperture. For example, we can construct the following URA sequence of length 6 . 

Example. PCSS with parameter (M, A, L) = (4, 2, 6 ): 

.Si = {1,-1,-1,-1,-1,-1}, S2 = (1,1,1,1,-1,-1}, S3 = {-1,1,-1,1,-1,-1}, .54 = {-1,-1, 1,1,-1,1}. 

Then, the sequences may be mapped onto a two-dimensional lattice, following procedures similar to Equations (fT4l i and ( fTbl l. 
Now a natural question that arises is: for a given alphabet, what are the possible lengths for which there exists a PCSS? and 
how to design them? This will be addressed in the remaining sections. 
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A natural way to construct PCSS is to synthesize them from existing designs. Some synthesis methods have been provided 
for binary PCSS in HO), and they could be easily extended to the non-binary case. In the following two sections, we propose 
some different synthesis methods. 

At the end of this subsection, there are two remarks worth mentioning. First, the concept of PCSS is not new. It was 
once referred to as “periodic complementary sequences” or “periodic complementary binary sequences” HOl . To the best of 
our knowledge, prior works mainly focused on the binary case. One possible reason is its intimate relationship with cyclic 
difference sets. Second, complementary sequence sets are subclasses of PCSS due to the following fact: 

Ap{v) = A^{v) + A^{v — L) VuGZ, 0<'!;<L, (17) 

where S' is a sequence of length L, Ap(-), Aa(-) respectively denote periodic, aperiodic autocorrelations. 

B. Synthesis Methods from the Chinese Remainder Theorem 

1) PCAS synthesized from PCSS and perfect sequence: A sequence is called a “perfect sequence” if its periodic autocorre¬ 
lation is a periodic extension of the discrete delta function. Consider a PCSS {Sm}m=i of length s, and a perfect sequence 
S of length t. We can then construct a PCAS {Cm}m=i of s x f (or similarly t x s): 

Cm[iM] = Sm[i]S[j], m = I,--- ,M, G Z. (18) 

Proof The periodic autocorrelation of Cm satisfies 

A^^{vi,V2) = A^''^{vi)A^{v2). 

Thus, for any vi,V 2 € Z,{vijV 2 ) ^ (0,0), 

M / M \ 

A^'^{vi,V2) = f Y I = 0. 

m—1 \m—l / 

□ 

2) PCSS synthesized from PCSS and perfect sequence: Consider a PCSS {Sm}m=i of length s, and a perfect sequence S 
of length t. Also assume that s and t are co-prime. We can then construct a PCSS {Sm}m=i of length st: 

= Cm[i mod mod f], i G Z, (19) 

where {Cm} is given in Subsection IIV-B 1 1 

Proof Equation ( fT9] l provides a one-to-one mapping between a sequence and an array, guaranteed by the Chinese Remainder 
Theorem. The mapping is linear so that the autocorrelation function is preserved, i.e. 

A^'^{v) = A^'^{v mod s,v mod t), 

and thus the sequence set {Sm}m=i complementary. □ 

3) PCSS/PCAS from two PCSS with co-prime lengths: Consider a PCSS {*S'mi}m^=i of length s, and another PCSS 
{Tm 2 }ml=i of length t. We can then construct a PCAS {C(^mi,m 2 )} of size s x t: 

C(mi,7n2)[hj] = Smi\i]Tm2[j], TOl = I,--- ,Mi,TO2 = 1, ' ' ' , M 2 , i,j G Z. (20) 

Further, if s and t are co-prime, we can construct a PCSS of length st. 

Proof For a given 1 < m 2 < M 2 , 

Ml 

Y. A^^"‘T-’^^'>(vi,V2 

mi—1 

Thus, for vi,V 2 G Z, (ui, U 2 ) f (0, 0), 

M 2 Ml 

Y^ Y^ 

7712 = 1 mi = l 

If s and t are co-prime, a PCSS could be designed using the mapping given in Equation ( fTOl l. □ 

Remark 8. This result is stronger than that given in 4401/ (Theorem 6), since it does not require the number of sequences to 
be relatively prime. 




s ■ (112) r;i = 0 


0 


otherwise 


M 2 / Ml \ 

m2 —I \mi—l / 
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4) PCAS constructed from another PC AS of a different size: Assume that we have a PC AS of size s x t synthesized from 
PCSS using the method in Subsection IIV-B3I Suppose that gcd(s,f) 7^ 1, but s = S1S2 for some 

Si 7^ 1 and S2 7^ 1 where gcd(si,S2) = gcd(s2,f) = 1. A PCAS of size si x S 2 t could be designed by first mapping the 
PCSS {Smi}ml=i ^ PCAS of size si x S2 in a way similar to Equation (fTSl l. then constructing a three-dimensional PCAS 
of size Si X S2 X f in a way similar to Equation (fTsT l. and finally mapping the latter two dimensions to a single dimension in 
a way similar to Equation ( fT9] l. resulting in a PCAS of size si x S 2 t. 


C. Synthesis via Unitary Matrices 

Theorem 5. For any positive integer s, there exists at least one PCSS with design parameters {M, N, L) = (p„, pi ■ ■ ■ pm s), 
where Pi < ■ ■ ■ < Pn tire all the distinct prime divisors of s. 

For any positive integers Si, • • • , s^, there exists at least one PCAS of size Si x • • • X s^ with design parameters (M, N, L) = 
{PmPi • • 'Pm st), where pi < ■ ■ ■ < p„ are all the distinct prime divisors 0/ Si • • • s^. 


Proof We prove the first part constructively. Using Equation dnii, we observe that it suffices to construct an aperiodic 
complementary array set. Without loss of generality, assume that s = pf x p® • • • x p‘f' is a prime factorization of s, 
where qj > 1 and pj’s are distinct for j = 1, • • • ,n. Let 5"^°^ = {<S'm^}m=i b® ^ set of pi sequences each of which contains 
a single point 1, i.e. has design parameters {M,N,L) = (pi, 1,1). 

In the first iteration, we apply Theorem[3to 5^°^ with U equal to the Eourier matrix Fp^ while satisfying Remark|3]conditions 
to obtain a (one-dimensional) complementary array set with {M,N,L) = (pi,lcm(pi, l),pi x 1) = (pi,Pi,Pi). Applying 
Theorem 3 a second time, we obtain a complementary array set with {M,N,L) = (pi, lcm(pi,pi),pi x pi) = (pi,Pi,p|). 
After applying Theorem 3 to 5'^°^ ?i — 1 times, we obtain the complementary array set with (M, N, L) = (pi,Pi,Pi^~ ). 

In the second iteration, we first apply Theorem 0] to with M = p 2 while satisfying Remark |5] conditions. The 
resulting complementary array set has parameters {M^N,L) = (p2,lcm(pi,p2),pi x p‘f~^) = (j)2,PiP2,Pi^)\ then we 
apply Theorem 3 q 2 — I times with U being the Eourier matrix Fp^ to create the complementary array set with 
{M,N,L) = {,P2,PiP2,pfpf~^)- 

By recursive construction as above, after ruth iteration we obtain the complementary array set with (M, N, L) = 
(Pn,Pi • • ■Pn,pTp '2 ' ■ 'Plr~^)- Finally, applying Theorem [3] with U equal to the Eourier matrix Fp^ an extra time to 
we obtain a complementary array set with (M, iV, L) = {pn,PiP2 ■ ■ 'PmP’1 P'2 '' 'P’n)- 

The proof of the second part is similar. 

□ 


Remark 9. Theorem^gives the construction for PCAS of an arbitrary size, where the alphabet is determined by the product 
of the distinct prime divisors of the size, and the order is the largest prime divisor. From the proof of Theorem |5] the result 
also holds for aperiodic autocorrelations. 

A natural question that arises is how tight the result in Theorem |5] is. Specifically, is there any solution whose order N is 
less than pn? This is clearly not the case for Golay complementary sequence pair, where the size s is a power of 2. Although 
we were not able to answer this question in general, we were able to prove the following results. 

Theorem 6. For any prime number p, a p-regular set is defined to be a set of p distinct unimodular complex numbers that 
form the vertices of a uniform polygon in the complex plane. Let N = p^f ■ ■ ■ pj^" ,rj > 1 be a positive integer with distinct 
prime divisors Pj,j = 1, • • • ,n. Consider M variables Xi, - ■ ■ , Xm that take values in the set of Nth root of unity. Suppose 
tttm = 0 . 

1) If n <2, the set {xm}m=i ttctn be written as the unions of p^-regular configurations, i.e. 

{xi,---,xm}= U (21) 

{(fe,i)|cfc>0, fe=l,2, ,Ck} 

2) M can be written as 


M = cipi-b C 2 P 2 , ci,C2eNU{0}. (22) 

Remark 10. Consider an aperiodic complementary array set {Sm}m=i with design parameters {M, N, L). Then we have 
Sm[0]Sm[L — 1] = 0. Assume that the alphabet is N-phase, where N = pf (n = \) or N = pfplf with pi and p 2 
distinct primes (n = 2). Applying Theorem^ Equation ( 1221 ) implies that (1) M >pi if n = 1; (2) M > min{pi,p 2 } if n = 2. 
Due to similar reasons, M = 7 in Design [7| is tight whenever N is a power of 7. 
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V. Simulation Results 

We have performed computer simulations to demonstrate a multi-channel CAI system and a classical URA-based one. 

The multi-channel CAI system that we select comes from Design |2] Admittedly, in practice we only need four pairs of 
aperture arrays with { — 1, l}-alphabet. But the coded images contain negative entries, which are not straightforward to illustrate 
by simulation (we used Matlab software). We thus provide an alternative approach which relies on the following lemma. 

Lemma 2. Suppose that {{Cm, Dm)}m=i complementary bank with alphabet 21 = {—1,1} and it satisfies 'Y2m=i = 0. 
Then {(Cm, Dm)}m=i ^ complementary bank, where Cm = (Cm + •^)/2, m = 1, • • • , M. Here, 0 and J are respectively 
the array of zeros and the array of ones, whose supports are the same as Dm- 

Proof The proof follows immediately from 

Y,Cm*Df^=Y,-{Cm+J)*Df^ = -Y,Cm*Dm + Z^J*[Y2Dm\ = ^ ^ ^23) 

m—1 m—1 m—1 \m—l / m—1 

□ 


The above result gives a general method to design a mask with simple closing/opening pinholes (the elements of C are 
either 0 or 1). In practice, the method may be of some interest on its own right, but we do not elaborate here. As a corollary 
of Lemma |2] it is easy to see that if {Cm}m=i is a complementary array set with alphabet 21 = {—1,1}, then 

+j),c'™^| y|Q(-c^ + 



is a complementary bank. 

Following Design |2] and the above result, we obtain the following 8 -channels CAI {{Cm, Dm)}m=i’ with a mask as 
shown in Figure |5] The coding arrays are 


{C'i[fc]}Lo = {1,1,0,1,1,0,1}, 

{C'4[fc]}Lo = {1,1,1,0,1,0,1}, 

{C'7[fc]}Lo = 10,0,1,0,1,0,0}, 


{C2[fc]}Lo = 10,1,1,0,0,1,1}, 
{C5[fc]}Lo = 10,0,1,0,0,1,0}, 
{C8[fc]}Lo = 10,0,0,1,0,1,0}. 


{C'3[fc]}Lo = (1,1,0,1,0,1,1}, 
{C'6[fc]}Lo = 11,0,0,1,1,0,0}, 


If we choose the element labeled 0 to be the origin, the decoding arrays are 


{i^iWlLo = 

{i^4[fc]}Lo = 

{D7[k]}l=o = 


( 1 , 1 , - 1 , 1 , 1 , - 1 , 1 }, {D^mi^o = 1 - 1 , - 1 , 1 , 1 , 1 , 1 , - 1 }, {D^mLo = ( 1 , - 1 , 1 , 1 , 1 ,- 1 , 1 }, 

{1,1,-1,1,1,1,-1}, {f?5[fc]}6^0 = 1-1,-1,1,-1,-1, 1,-1}, {Z?6 [A]}Lo = {1, 1,-1,-1,-1,-1,1}, 

1-1,1,-1,-1,-1,1,-!}, {i^8[/c]}Lo = 1-1,-1,1,-1,-1,-1,1}. 


Fig. [T 6 ] illustrates how a source object is coded and decoded in a multi-channel system. The source object is a 131 x 131 pixels 
“camera man” image and the distance between two pinholes is 60 pixels. The grey level is normalized to be in the range of 
[0,255]. 

The second simulation is for URA (shown in Fig. [TtI i. It uses the same source object and following aperture: 


0 1111 
0 0 110 
0 10 0 1 


(24) 


which is a 3 X 5 array that comes from the m-sequence of length 15 (See Section UlI-AI for details). In the implementation, 
we use the arrangement suggested in 0 , i.e. a 6 X 10 aperture composed of a periodic extension of the basic 3x5 patterns, 
and a 3 X 5 decoding array. 


Appendix A 
Proof of Theorem[T] 

The following lemmas are helpful to the proof of Theorem [T] 

Lemma 3. If complex numbers cti, • • • , a 4 are unimodular and satisfy X]t=i = 0, then they contain 2 opposite pairs. 
Proof of Lemma |i| 

Let y = Qfi -I- a2, z = (—0:3) + (-0:4), then y = z. Because |ai| = \a2\ = 1 , y is on the bisector of ai and a2- Similarly, 
z is on the bisector of —0:3 and —a^. Since y = z. We have ol\ = —0:3 or a\ = -04. 

Lemma 4. If roots of unity ai,a 2 satisfy \ai — 2a2\ = 1, then ai = a 2 - 
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(a) The upper-left image, “cameraman”, is the source image. From the upper-middle to the bottom-right, the images are the coded images from 
apertures Ci, • • • , Cs in each channel. 



(b) From the upper-left to the bottom-middle, the images are the decoded results from Di, • • • , Dg in each channel. The bottom-right image 
is the reconstmcted image, coming from the addition of the eight decoded results. 

Fig. 16. Demonstration of the encoding and decoding process of a multi-channel CAI system 
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(a) This figure shows the coded image of the URA system. The coded aperture produces a cyclic version 
of the basic aperture pattern. 



(b) This figure shows the decoded image. The source image is reconstructed in the center area. 
Fig. 17. Demonstration of the encoding and decoding process of a URA-based CAI system 
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Proof of Lemma 0 

The identity 1 = (ai — 2 a 2 ){ai — 2 d 2 ) = 1 — 2 (aid 2 + dia 2 ) + 4 implies aid 2 = 1, i.e. cti = a 2 - 
Proof of Theorem [7J 

Assume that there exists a complementary array pair. Writing down (l2]i explicitly we obtain the following system of equations 


consisting of 9 equations and 14 variables {xk}%^Q U {?/fc}fc=o- 

xixz + x^XA + yiVd. + y&VA = 0 (25) 

ato^i + X'iX 2 + xsxe + xaxq + yoyi + ysm + 2/5^6 + yiVo = 0 (26) 

xiXA + yiTjA = 0 (27) 

X 2 XA + xix^ + 2/22M + yiTJb = 0 (28) 

a;o ^2 + XqXi + xax^, + + 2 / 02/2 + 2/62/T + 2/42/“3 + 2/52/“o = 0 (29) 

a;2^5 + 2/22/“5 = 0 (30) 

a;3^“5 + X 2 XQ + 2/32/“5 + 2/22/”6 = 0 (31) 

x^xz + X1X2 + x^XA + Xexo + yoys + 2/12/2 + 2/52/“4 + 2/62/”o = 0 (32) 

a;3^6 + 2/32/“6 = 0 (33) 

We only need to prove that the above system of equations have no solution on the unit circle. Assume without loss of generality 
that xi = 2/1 = 1- After simplifying Equations (l27l i. (l25l l. (l30l l. (l3^ we have 

2/4 = -XA, (34) 

2/3 = -( 3^3 + (^6 - ^4)a;4) = -X 3 - xqXa + 1, (35) 

2/5 = -X2X5y2, (36) 

2/6 = -xsxeys. (37) 

From Equation (l35T l. we obtain ys + X 3 + xqXa — 1 = 0. Due to Lemma |3 we have three cases to consider; 

Case A: 2/3 = l,a ;4 = —X 3 Xe', Case B: X 3 = 1 , 2/3 = —xaxq'. Case C; xq = XA,y 3 = —X 3 . 


A. Case A: y 3 = 1,xa = —X 3 Xe 

From Equations ( |34] |. ( |36] |. and (IJTI i. we eliminate variables X 4 , 2 / 3 , 2 / 4 j 2 / 51 2/6 in Equations ( |28] ) and OTT ) and obtain 

-X 2 X 3 XQ + *5 + X3XQy2 - X 2 X 3 y 2 = 0 (38) 

X 3 X 5 + X 2 Xe - X 2 X 5 y 2 - X3Xey2 = 0. (39) 

We may write Equations ( |38] | and ( |39] | as 

a;5^3(a;2 - 2 / 2 ) = a;6(l - a; 22 /“ 2 ) (40) 

0:5 {X3 - X2y2) =Xe {X3 2/“2 - X 2 ) • (41) 

If X 2 y 2 , Equation (l40l i gives 

1 - X2y2 

X 5 = X 3 - X 6 = -X 3 y 2 Xe. 

X 2 - 2/2 

If X 2 X 3 y 2 , Equation dTO gives 

X 3 V 2 — X 2 

X 5 = — - ^-^Xe = 2 / 2316 - 

X 3 - 3122/2 

So there are four cases to consider that further eliminate the variables. 

Case Al: X 5 = - 3 : 32 / 2 X 6 , X 5 = fzxe 

Clearly, X 3 = —1. Eliminating variables in Equations (l26l l and (l3^ we obtain 

(xo - x “2 - x “2 - xeyo) + ( 2/2 + xex'o + 2/0 + m) = 0 
-(xo - X 2 - X 2 - Xeyo) + ( 2/2 + xexb + 2/0 + 2 / 2 ) = 0 


which is equivalent to 


Xo - X2 - X2 - Xeyo = 0 
2/2 + X 6 X“o + 2/0 + 2/2 = 0 . 


(42) 

(43) 




21 


We further obtain 

yo = X(i{xo - 2x2) (44) 

y2 = xe{x2 - xq). (45) 

Equation (l44l) gives jafo — 2 x 2 1 = 1, which further implies that xq = X 2 - This is a contradiction to Equation (l45l l. 

Case A2: X 2 = y 2 ,X 5 = y 2 Xe 

We only need to check the validity of Equations (|2^ . (|2^ . and (l32] |. We write them in terms of 5 variables xq, X 2 , X 3 , xe, yo: 

Xq + X3X2 + X2- xsXqXo + t/o + ^2 + ;c 2 a ;3 + X^X^yQ = 0 ( 46 ) 

xoX2 + X2X6X0 + yoX 2 - X3X6 + X3X6 - X 2 Xeyo = 0 

Xq + X3X2 - X2 + X3X6X0 + yoX 3 + X2X3 - X2 - xeyo = 0 . ( 47 ) 

In fact, a contradiction can be obtained from Equations (|46] | and (l47l) . Taking the sum and difference of the two equations, we 
obtain: 


(x3 + l)yo + (x3 - l)xeyo = -4x2X3 - 2xo (48) 

(X3 + l)x6t7o - (X3 - l)yo = -4X2 + 2 x3X6X“o. (49) 

If we have a valid solution (x2, X3, Xq, yo; it is easy to see that {^~^X 2 ,X 3 ,^XQ,^yQ,^‘^Xe) is also a valid solution for 

any unimodular complex number Therefore, we only need to consider the case xe = 1. Replacing Xe = 1 into Equation 

gUl, multiplying the equation with —X3, and then taking the conjugate, we obtain 


- (X 3 + 1)2/0 - { X3 - 1)2/0 = 4X2X3 - 2 X 0 - 


Adding Equations (|48] | and (ISOl l gives 


a^o = (2:2 - a:2)x3. 


Erom the identity 1 = |xo| = |x2 — X2IIX3I = |x2 — X2I, X2 is in the form of 

^^1 + 2 ^ 2 , (51,^2 G {1,-1}- 

Combining Equations (l52l i and (BTl i gives xq = 162 X 3 . Eurthermore, Equation (|48]) is simplified to be 

(X3 + 1)2/0 + (2:3 - l) 2 fo = -2V3S1X3. 


Equation ( 1531 ) implies that 

2-^3 = |(x3 + 1)2/0 + (2:3 - l)2/b| < |a;3 + 1| + |a;3 - 1| < 2V2, 


(50) 

(51) 


(52) 

(53) 


which is a contradiction. 

Case A3: X5 = -X32/2X6, X2 = X32/2 

Eirst, we rewrite Equations (l26l l and ( l32l i in terms of xq , X3 , xe , 2/0: 2/2 ^ 


2:0 + 21/2 - 2x32/2 - X3XeXo + 2/0 + xsxeyo = 0 (54) 

X 0 X 3 + 2x32/“2 + 22/2 + X 6 x“o + 2/0 - X3Xeyo = 0. (55) 

If we have a valid solution ( 2 / 2 , 2 : 3 , xq, 2/Oj aie), it is easy to see that {^~^y 2 ,X 3 ,^XQ,^yo,^^xe) is also a valid solution for 
any unimodular complex number Therefore, we only need to consider the case xo = 1. By computing Equation ( |5^ + X 3 - 
Equation (fSSl) and Equation (l5^ - X 3 - Equation (l55l) we obtain 

(X3 + 1 ) 2/0 + ( 2:3 - l)2/“o = -42/2 - 2xo (56) 

-(x3 + l)2/“o + (x3 - 1 ) 2/0 = -4x32/2 - 2 x3X“o. (57) 

Multiplying Equation (l57l i by X 3 , and then taking the conjugate, we obtain 

- ( 2:3 + 1 ) 2/0 - ( 2:3 - 1 ) 2/0 = -41/2 - 2xo. (58) 

Adding Equations ( 1561 ) and ( 1581 ) we obtain 

2:0 = -( 2/2 + 2 / 2 )- 


2/2 = 2'tl + “^"2, 


61,62 G {1,-1}. 


Thus, 2/2 is in the form of 
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Furthermore, 

Xo = -(5i 

(X 3 + l)yo + {X 3 - l)yo = -2V3S2, 

which implies that 2\/3 < 2\/2. 

Case A 4 : X2 = y2,X2 = 2 : 32/2 

Similar to case Al, Equations (l26l l and (l3^ imply 

-{2x2 + xo) = yo (59) 

X5Xe - xeixo + X 2 ) = 0- (60) 

Equation ( l59b implies that X 2 = —Xq, which is a contradiction to Equation ( l60b . 


( 61 ) 

( 62 ) 

gives X 5 = —X 42 / 2 . We therefore have the 


B. Case B: X3 = l,?/3 = —x^Xq 

Erom Equation ( l 28 l l and OTT i. we obtain 

(®5 - y2X4)il - X 2 y 2 ) = 0 
{X 2 X 3 + 2 / 2 ^ 4 ) (1 + 2 : 4 ^ 52 / 2 ) = 0 . 

If 2/2 X4X3, Equation ( 1 ^ gives X2 = 2/2- If X2 ^ —X6X42/2, Equation ^2 

following four cases to discuss. 

Case Bl: 2/2 = 2:4X5, X2 = —X6X42/2 = —X3X3 

Eirst, we rewrite Equations (| 29 ] | and ( | 32 ] | in terms of xq, X4, X5, xe, 2/0^ 

2 x 6 - X0X5X6 + 2 x 4 + X4X52/0 + X4X5X62/0 + X0X5 = 0 ( 63 ) 

2x4x5 + Xo - 2x5X6 - X4X62/0 + X42/0 + X0X6 = 0 . ( 64 ) 

If we have a valid solution (xq, 2/Oj X4, X6, X5), it is easy to see that (^Xq, ^2/0; C^X4, ^^xe, ^^X5) is also a valid solution for 
any unimodular complex number Therefore, we only need to consider the case xe = 1. Taking X6 = 1 into Equations (|6^ 
and dMl i. multiplying Equation ( l 64 l i by —X5, and taking its conjugate, we obtain 

2 - X0X5 + 2 x 4 + X4X52/0 + X 4 X 52 /b + xqXs = 0 ( 65 ) 

2 - X0X5 - 2 x 4 - X4X52/0 + X 4 X 52 /b - X0X5 = 0. (66) 

Adding Equations (l 65 l l and (l6^ gives 

4 = 2x0X5 - 2x4X52/0. ( 67 ) 

Because 1 2x0X5] = 12x4X52/0 1 = 2 , the only possibility is 

X 0 X 5 = l,X4X52/b =- 1 - ( 68 ) 

Taking Equation (l68l l into (l 65 T l we obtain 2x4 = 0 , which is a contradiction. 

Case B 2 : X2 = 2/2,X2 = —X6X42/2 

Clearly, xe = —X4. Because 2/3 = —X4X6 = 1 , this case is covered by Case A. 

Case B 3 : 2/2 = X4X5,X5 = —X42/2 

Clearly, 2/2 = —X4X5. Because 2/2 = X4X5 = —2/2,This case is not possible. 

Case B 4 : 0:2 = 2/2, X5 = —X4272 = —X4X2 

Eirst, we rewrite Equations (l 26 l l and ( l 29 l l in terms of Xo, X2, X4, xe, 2/0^ 

2 x 2 — 2X2X4X6 + Xo + 2/0 + X4X0 — X 427 o = 0 ( 69 ) 

2 x 4 + 2 x 6 + X0X2 + X22/0 — X2X4X0 + X 2 X 427 o = 0 . ( 70 ) 

If we have a valid solution (x2, xq, 2/Oj X4, xe), it is easy to see that (^“^X2, ^xq, ^2/0i C^X4, ^^X6) is also a valid solution for 

any unimodular complex number Therefore, we only need to consider the case X4 = 1 . Taking the conjugate of Equation 

(| 69 ] |, multiplying Equation (iTOl i by X2, and letting X4 = 1 , we have 

2X2 - 2X2X6 + Xq + 2 /“o + Xo - 2/0 = 0 ( 71 ) 

2 x 2 + 2 x 2 X 6 + Xo + 2/0 - Xo + 2 /b = 0 . ( 72 ) 

Adding Equations ItTI i and (l 72 l i. we obtain 

4x2 + 2 xo + 22/0 = 0 . ( 73 ) 

Because |4x2| = 4 , | 2 xo| = | 227 o| = 2 , the only possibility is 

Xo = -X 2 ,yo = -X 2 . 

Applying Equation (| 7 ^ to (TtTT i gives 2x2X6 = 0 , which is a contradiction. 


( 74 ) 
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C. Case C: Xq = X4,y3 = —X3 

From Equations (l 28 l l and OTl i. we obtain 

{X5 - X 4 y 2 ){i - X2y2) = 0 ( 75 ) 

{X3X3 + X4y2){x2 + 2 / 2 ) = 0. ( 76 ) 

If 2/2 ^ X2, Equation (l 75 l l gives ?/2 = X4X3. If 2/2 ^ —X2, Equation 17 ^ gives ?/2 = —X3X4X5. So there are four cases to 
consider. 

Case Cl: 2/2 = X2, 2/2 = —X3X4X3 

Eirst, we rewrite Equations (| 29 ] | and (l 32 l i in terms of xq, X 2 , X 3 , X 4 , t/o: 

2 x 4 + 2x3x4 + ^20:0 + X2yo - X2X3X4X0 + X2X3X4yo = 0 ( 77 ) 

2 x 2 — 2.X2X3 + XcXq — X 3 yo + X4X0 + X 4 yo = 0. ( 78 ) 

If we have a valid solution {x2,X3,xo,yo, X4), it is easy to see that xa, ^ 2 / 0 ; ^^^ 4 ) is also a valid solution for any 

unimodular complex number Therefore, we only need to consider the case X 4 = 1. Taking the conjugate of Equation ( iTSl l, 
multiplying it by X2, and letting X4 = 1, Equations (ITTT i and ( fTHI l give 


2 + 2 x 3 + X 2 X 0 + X 22/0 - X 2 X 3 X 0 + X 2 X 32/0 = 0 

(79) 

2 - 2x3 + X 2 X 0 + X 22/0 + X 2 X 3 X 0 - X2X3yo = 0. 

(80) 

Adding Equations (|7^ and dSOll gives 

4 + 2x2 (xo + 2 / 0 ) = 0. 

(81) 

Thus |xo + 2 / 0 I =2, the only possibility is 

Xo = 2/0 = -X2- 

(82) 

Applying Equation (l82li to (l79ll, we obtain 2 x 3 = 0, which is a contradiction. 

Case C2: 2/2 = X 4 X 5 , 2/2 = —X 2 

Eirst, we rewrite Equations (l26ll and (|29]| in terms of Xq, X 2 , X 3 , X 4 , 2 / 0 : 


2 x 3 X 2 - 2 x 2 + Xo + 2/0 + X 4 X 0 - X42/b = 0 

(83) 

2 X 4 + 2 X 4 X 3 + X 0 X 2 — 2 / 0 X 2 — X 2 X 4 X 0 — X2X42/b = 0. 

(84) 

If we have a valid solution (x 2 , X 3 , xo, 2 / 0 , X 4 ), it is easy to see that (^“^X 2 , X 3 , ^xo, ^ 2 / 0 , C^X 4 ) is also a valid solution for any 
unimodular complex number Therefore, we only need to consider the case X 4 = 1. Taking the conjugate of Equation (l84li. 
multiplying it by X 2 , and letting X 4 = 1, Equations (l83Tl and (l84li give 

2 x 3 X 2 - 2 X 2 + Xo + 2/0 + - 2/0 = 0 

(85) 

2 x 3 x 2 + 2 x 2 - Xo - 2/0 + Xo - 2/0 = 0. 

(86) 

Subtracting Equations dSbll and dSSll. we obtain 


4X2 - 2(xo + 2 / 0 ) = 0. 

(87) 

Thus |xo + 2 / 0 I =2, and the only possibility is 

Xo = 2/0 = X“2. 

(88) 


Applying Equation (l 88 l l to (l85l l. we obtain 2 x 3 X 2 = 0, which is a contradiction. 

Case C3: 2/2 = X 4 X 5 , 2/2 = — X3X4X5 This case is covered by Case A, because X3 = — 1 , 2/3 = — X3 = 1. 

Case C4: 2/2 = a^ 2 , 2/2 = —X 2 This case is clearly not possible. 


Appendix B 
Proof of Theorem[2] 

The proof follows a similar procedure to that of the proof of Theorem [T] The only difference is that the modulus of xq , 2/0 
are changed from one to zero. Related changes in the proof in Appendix are listed below. 

Case Al: Equation (l42T i gives 2 x 2 = 0, which is a contradiction. 

Case A2: Equation (l48l l gives 4 x 2 X 3 = 0, which is a contradiction. 

Case A3: Equation (l5^ gives 42/2 = 0, which is a contradiction. 

Case A4: Equation (l5^ gives 2 x 2 = 0, which is a contradiction. 

Case Bl: Equation dCTl i gives 4 = 0, which is a contradiction. 

Case B4: Equation (l73l l gives 4 x 2 = 0, which is a contradiction. 

Case Cl: Equation dSTl i gives 4 = 0, which is a contradiction. 

Case C2: Equation dSTl i gives 4 x 2 = 0, which is a contradiction. 
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Appendix C 
Proof of Theorem[6] 

Let Z[A] denote the polynomial ring over Z, and ^n{^) denote the iVth cyclotomic polynomial. Let ^ = exp(*27r/iV), 
Un = I j = 0, • • • , TV — 1} be the group of TVth roots of unity endowed with multiplication. For any r] G Un, let \ri\ 
denote the order of p in the cyclic group Un. Because Ylim=i there exists a polynomial F{\) = ^ 

such that fj > 0 and F(^) = 0 . 

We prove Theorem |6] using a sequence of lemmas. 

Lemma 5. Let be distinct prime numbers and integers > 0, k = 1,2. Then 

<l>p,(At’i^"'t’?) = n<i>p^Xp.(A), (89) 

i^O 

(90) 

$p,(Af"") = $p^.(A). (91) 

Similar results hold if pi and r\ are respectively replaced with p 2 and r 2 in the above equations. 

Proof of Lemma |5] 

Since both sides are monic and have degree (pi —= (pi —it suffices to show that every zero of <i> ri (A) 

r — 1 7’— 1 -1 

is a zero of $pj(Afi ). If p is a zero of $ r-i *2 (A), then |p| = pf^p^f, which implies |pt>i | = |p|/gcd(|p|,p^^“ ) =pi. 

Therefore, p is also a zero of $pj^(Afi ). 

The proof of Equations ( |90] | and (IgTT i is similar. 

Lemma 6 . If n = 1, i.e. N = p'f, then there exists a polynomial A(A) = ^ G Z[A], aj > 0, such that 

Proof of Lemma [d] 

First, <I’Ar(A) divides F{X), because F{X) annihilates ^ and $Ar(A) is an irreducible and monic polynomial in the ring Z[A]. 
Besides this. Equation dMI ) gives (A^/fi) = <i)jv(A). Therefore, there exists a polynomial Afc(A) = ^ 

such that F(A) = 

Second, because deg(E') < TV, we have deg(A) = deg(E') — deg(<hAr(A)) < Pi^~^ = TV/pi. We note that fj = aj,j = 
0, • • • ,TV/pi — 1 and that fj > 0. Therefore, the coefficients of A(A) are nonnegative. 

Lemma 7. If n = 2, i.e. TV = p^^plf, then there exist polynomials A],{X) G Z[A], fe = 1, 2 such that 

ii(A)$p,(A^) + i 2 (A)$p,(A^) = $jv(A). (92) 

Proof of Lemma [71 

First, Equation dS^ and its similar result (by replacingpi, ri with P2,72) imply that gcd($pj [X^/p^), ‘I’pj (A^/^^)) = <i>Ar(A). 
Second, consider two polynomials Tt^ (A) = 1 + A + • • • + k = 1 , 2 , where ti > t2 and gcd(Ti, ^2) = 1 - We apply 

Euclidean division to Ti, ^2 to obtain ti = t2q + b, 0 < b < t2. It is easy to observe that (A) = Tt^ (A) X*^~A2 +Ti,(X). 

If we continuously apply Euclidean division, we will find polynomials Ak{X) G Z[A], fc = 1 , 2 such that 

ii iX)Tt, (A) + i 2 iX)Tt, (A) = 1 (93) 

Replacing ti,t 2 , and A respectively by pi,p 2 , and A^/^t’iP 2 ) jjj Equation (i93T l, multiplying both sides by $ 7 v(A), and using 
Equation ( i90l l and its similar result, we obtain Equation (i92] i. 

Lemma 8 . If n = 2, i.e. TV = Pi^P 2 ^, then there exist polynomials Afe(A), deg(Afe) < N/pk — l,k = 1,2 such that 
FW = El=iM>^)HkiX), where HkW = $p,(A^/P'=), fc = 1,2. 

Proof of Lemma 0 

Clearly, <hAr(A) divides F{X) due to the reason mentioned before. Therefore, Lemma [7] implies that there exist polynomials 
Ak{X) G Z[A], fc =1,2 such that F{X) = X]fc=i -^k{X)IIk{X) holds. 

It is easy to see that X‘^Hk{X) can be written as X‘^Hk{X) = (A^ — 1)(5(A) + X‘^°Hk{X) for some Q{X) G Z[A],0 < do < 
N/pk — 1. Thus, there exist polynomials Ak{X), deg(Afc) < N/pk — 1 , k = 1,2, and IL(A) such that 

2 

F(A) = ^ Ak{X)Hk{X) + W{X){X^ - 1). 

k^l 

Since deg(F) < N/pk — 1, we obtain VF(A) = 0. 
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Lemma 9. The coefficients of Ak{\), k = 1,2 in Lemma \^can be made nonnegative. 

Proof of Lemma [P] 

Let 


PI 


--1 


P2 


--1 


— < [oio, • • • ai(7v/pi-i), 020, • • • , 02(Ar/p2-i)] I ^ ^ a2j2H2{X). 


21=0 

_N_i 


22=0 

_N_1 




( 2 ) 


N 


[oio, • • • ai( 2 v/pi_i), 020 , • • • , a 2 (Ar/p 2 -i)] | F{X) — ^ aij^Hi{X) + ^ a 2 j 2 -ff 2 (A), 02^2 > 0 ,^2 — 0, • • • , — — 1. 


21=0 


22=0 


For a fixed integer 0 < j < N/pi — 1, consider the set {j + kN/pi, fc = 0, • • • ,pi — 1} (mod N). For each k = 0, - ■ ■ ,pi — l, 
we apply Euclidean division to j + kN/pi and N/p 2 , and obtain integers 0 < pt < N/p 2 — 1,0 < hk < P 2 — ^ such that 
j + kN/pi = Pk + h}.N/p 2 . Because of the identity 

P 2-1 Pi-i 

IJ {fffc + {hk+T)N/p 2 ,k = 0, •• • ,pi - 1} = IJ {pfc + {hk + T)N/p 2 ,T = O,-- • ,p 2 - 1} (mod N), 

r =0 fc =0 

within the set Tip we can always decrease ai^j+rN/p^imoAN/pi)),''' = 0, • • • ,p 2 — 1 by one, while increasing a 2 gf,,k = 

( 2 ) 

0 , • • • ,pi — 1 by one. We conclude that the subset SL of Sir is not empty. 

To finish the proof, it suffices to show that within Dp , there exist an element with oy^ > 0 for all ji = 0, • • • , N/pi — 1. 
If this is not true, then there exists p, < 0 such that 


P = 


[aio,---ai(^,/p,,„i),a2o,-- - ,a2(N‘/p2-i)]S2)F 


( 2 ) 


I min {oio, • • • , Oi(Ar/pi_i)} | < 0 . 


(94) 


Suppose that aij = p. For each k = 0, ••• ,pi — 1, consider the nonnegative coefficient of the item in F{X): 

24i(A)iFi(A) contributes a negative value aij to it, and thus A 2 {X)H 2 {X) contributes a positive value. In other words, there 
exist integers 0 < pk < N/p 2 — 1,0 < hk < P 2 — 1 such that j + kN/pi = pk + hkN/p 2 and that a 2 gk > 0. It is clear 
that pk,k = 0, - ■ ■ ,pi — 1 are distinct values. By similar reasoning as before, we can increase ai(j+rAf/p 2 (mod M/pi)), = 

0 , • • • ,p 2 — 1 by one, while decreasing a 2 g,, ,k = 0, - ■ ■ ,pi — 1 by one, in order to get another element in Dp . Thus, we can 

increase | min {uio, • • • , } j, contradicting the definition of p in dH- 

Proof of Theorem | 6 ] ; 

Combining Lemmas | 6 ] to |9] we conclude that F{X) can be written as: 

-5L_1 

2 Pk 

pw = E ('^); where ^ ^ ^ ^kj ^ 0 : 

k^l j^O 

Hk{X) = 4>p,(A^/P'“) = 1 + A^ +A"^ + 

which is equivalent to Equation (| 2 TI) . Equation (l 22 l i then immediately follows from Equation (| 2 TI) . 
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